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We compare parton densities for heavy quarks. 



Reactions with incoming heavy (c,b) quarks are of- 
ten calculated with heavy quark densities just like 
those with incoming light mass (u,d,s) quarks are cal- 
culated with light quark densities. The heavy quark 
densities are derived within the framework of the so- 
called zero-mass variable flavor number scheme (ZM- 
VFNS). In this scheme these quarks are described by 
massless densities which are zero below a specific mass 
scale /x. The latter depends on m c or wib- Let us 
call this scale the matching point. Below it there are 
rif massless quarks described by rif massless densities. 
Above it there arc n.f + 1 massless quarks described by 
nf + 1 massless densities. The latter densities are used 
to calculate processes with a hard scale M 3> m c , m^. 
For example in the production of single top quarks via 
the weak process qi + b — > qj + t, where g,-, qj are 
light mass quarks in the proton/antiproton, one can 
argue that M = vtit should be chosen as the large scale 
and mb can be neglected. Hence the incoming bottom 
quark can be described by a massless bottom quark 
density. 

The generation of these densities starts from the so- 
lution of the evolution equations for n / massless quarks 
below the matching point. At and above this point 
one solves the evolution equations for nf + 1 mass- 
less quarks. However in contrast to the parameteri- 
zation of the ^-dependences of the light quarks and 
gluon at the initial starting scale, the x dependence of 
the heavy quark density at the matching point is fixed. 
In perturbative QCD it is defined by convolutions of 
the densities for the n / quarks and the gluon with spe- 
cific operator matrix elements (OME's), which are now 
know up to 0(a 2 s ) [ 0. These matching conditions 
determine both the ZM-VFNS density and the other 
light-mass quark and gluon densities at the matching 
points. Then the evolution equations determine the 
new densities at larger scales. The momentum sum 
rule is satisfied for the nt + 1 quark densities together 
with the corresponding gluon density. 
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Figure 1. The charm quark density :ecnnlo(5, x, fi 2 ) 
in the range 1CT 5 < x < 1 for \± 2 = 20.25, 25, 30, 40 
and 100 in units of (GcV/c 2 ) 2 . 

Parton density sets contain densities for charm and 
bottom quarks, which generally directly follow this ap- 
proach or some modification of it. The latest CTEQ 
densities [ ^) use 0(a s ) matching conditions. The x 
dependencies of the heavy c and b-quark densities are 
zero at the matching points. The MRST densities [ || 
have more complicated matching conditions designed 
so that the derivatives of the deep inelastic structure 
functions F2 and Fl with regard to Q 2 are continu- 
ous at the matching points. Recently we have pro- 
vided another set of ZM-VFNS densities [ f|, which 
are based on extending the GRV98 three-flavor den- 
sities in [ P to four and five-flavor sets. GRV give 
the formulae for their LO and NLO three flavor den- 
sities at very small scales. They never produced a c- 
quark density but advocated that charm quarks should 
only exist in the final state of production reactions, 
which should be calculated from NLO QCD with mas- 
sive quarks as in [ 0] . We have evolved their LO and 
NLO densities across the matching point /1 = m c with 
0(a 2 s ) matching conditions to provide LO and NLO 
four-flavor densities containing massless c-quark den- 
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Figure 2. Same as Fig.l for the NLO results from 

CTOfi e .„+ 1 



Fig. 3 
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Figure 4. The bottom quark density o;6nnlo(5, x > 
in the range 1CT 5 < x < 1 for fj? = 20.25, 25, 30, 40 
and 100 in units of (GcV/c 2 ) 2 . 

In the CS set they start at fi 2 = ml = 20.25 GeV 2 . At 
this scale the charm densities in the CS, MRST98 (set 
1) and CTEQ5HQ sets are shown in Figs. 1,2, 3 respec- 
tively. Since the CS charm density starts off negative 
for small x at /i 2 = m 2 = 1.96 GeV 2 it evolves less than 
the corresponding CTEQ5HQ density. At larger n 2 all 
the CS curves in Fig.l are below those for CTEQ5HQ 
in Fig. 3 although the differences are small. In general 
the CS c-quark densities are more equal to those in the 
MRST (set 1) in Fig.2. 
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Figure 3. Same as Fig.l for the NLO results from 
CTEQ5HQ. 

sities. Then these LO and NLO densities were evolved 
between \i = m c and fi = with four-flavor LO 
and NLO splitting functions. At this new matching 
point the LO and NLO four-flavor densities were then 
convoluted with the 0(a 2 ) OME's to form five-flavor 
sets containing massless b-quarks. These LO and NLO 
densities were then evolved to higher scales with five- 
flavor LO and NLO splitting functions. Note that the 
0(a 2 s ) matching conditions should really be used with 
NNLO splitting functions to produce NNLO density 
sets. However the latter splitting functions are not yet 
available, so we make the approximation of replacing 
the NNLO splitting functions with NLO ones. 

In this short report we would like to compare the 
charm and bottom quark densities in the CS, MRS and 
CTEQ sets. We concentrate on the five-flavor densi- 
ties, which are more important for Tevatron physics. 
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Figure 5. Same as Fig. 4 for the NLO results from 
MRST98 set 1. 

At the matching point /i 2 = 20.25 GeV 2 the b-quark 
density also starts off negative at small x as can be seen 
in Fig. 4, which is a consequence of the explicit form of 
the OME's in [ |. At 0(a 2 ) the OME's have non- 
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logarithmic terms which do not vanish at the match- 
ing point and yield a finite function in x, which is the 
boundary value for the evolution of the b-quark den- 
sity. This negative start slows down the evolution of 
the b-quark density at small x as the scale /x 2 increases. 
Hence the CS densities at small x in Fig. 4 are smaller 
than the MRST98 (set 1) densities in Fig. 5 and the 
CTEQ5HQ densities in Fig. 6 at the same values of /i 2 . 
The differences between the sets are still small, of the 
order of five percent at small x and large fi 2 . Hence 
it should not really matter which set is used to cal- 
culate cross sections for processes involving incoming 
b-quarks at the Tevatron. 
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to decide which set is the best one. We end by noting 
that all these densities are given in the MS scheme. 

REFERENCES 



2. 



M. Buza, Y. Matiounine, J. Smith, W.L. van Neer- 
ven, Eur. Phys. J. CI, 301 (1998). 
H.L. Lai, J. Huston, S. Kuhlmann, J. Morfin, F 
Olncss, J. Owens, J. Pumplin, W.K. Tung, licp- 
ph/9903282| . 



Roberts, W.J. Stirling and R. 
J. C4, 463 (1998). 



3. A.D. Martin, R.G. 
Thome, Eur. Phys 

4. A. Chuvakin, J. Smith, [hcp-ph/9911504] , to be pub- 
lished in Phys. Rev. D. 

5. M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C5, 
461 (1998). 

6. E. Laenen, S. Riemersma, J. Smith and W.L. van 
Neerven, Nucl. Phys. B392, 162 (1993); ibid. 229 
(1993); S. Riemersma, J. Smith and W.L. van 
Neerven, Phys. Lett. B347, 43 (1995); BW. Harris 
and J. Smith, Nucl. Phys. B452, 109 (1995). 



Figure 6. Same as Fig. 4 for the NLO results from 
CTEQ5HQ. 



We suspect that the differences between these re- 
sults for the heavy c and b-quark densities are pri- 
marily due to the different gluon densities in the three 
sets rather to than the effects of the different boundary 
conditions. This could be checked theoretically if both 
LO and NLO three-flavor sets were provided by MRST 
and CTEQ at small scales. Then we could rerun our 
programs to generate sets with 0(a 2 s ) boundary con- 
ditions. However these inputs are not available. We 
note that CS uses the GRV98 LO and NLO gluon den- 
sities, which are rather steep in x and generally larger 
than the latter sets at the same values of /i 2 . Since the 
discontinuous boundary conditions suppress the charm 
and bottom densities at small x, they enhance the the 
gluon densities in this same region (in order that the 
momentum sum rules are satisfied). Hence the GRV98 
three flavour gluon densities and the CS four and five 
flavor gluon densities are generally significantly larger 
than those in MRST98 (set 1) and CTEQ5HQ. Unfor- 
tunately experimental data are not yet precise enough 
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